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Abstract 

Using a Generic spin basis, we present a general formalism of one-loop radiative 
corrections to the spin correlations in the top quark pair production at the Next 
Linear Collider, and calculate the 0(a s ) QCD corrections under the soft gluon 
approximation. We find that: (a) in Off-diagonal basis, the 0(a s ) QCD corrections 
to efe + ( e^e + ) scattering process increase the differential cross sections of the 
dominant spin component M\ (t\tt) by ~ 30% and ~ (0.1% — 3%) depending on 
the scattering angle for = 400GeU and ITeV, respectively; (b) in Off-diagonal 
basis (Helicity basis), the dominant spin component makes up 99.8% (~ 53% ) of 
the total cross section at both tree and one- loop level for ^/s = 400GeU, and the 
Off- diagonal basis therefore remains to be the optimal spin basis after the inclusion 
of 0(a s ) QCD corrections. 



PACS numbers:12.38.Bx, 13.88.+e, 14.65.Ha 
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A special feature of the top quark is that due to its large mass, m t = 175 ± 6GeV[0], 
the lifetime of the top quark is very short. This has an important consequence that 
the top quark decays before it hadronizes@, and the spin information is preserved from 
production to decay. Thus we can expect the spin orientation of the top quark to be 
observable experimentally. Since there are significant angular correlations between the 
decay products of the top quark and the spin of the top quark, the spin-spin correlations 
in the top quark pair production can be extracted by forming angular correlations among 
the decay products of the top quark and the top anti-quark. Thus the spin correlations in 
the top quark pair production can be used as a good observable for testing the standard 
model (SM). 

In most papers on the spin correlations for top quark pair production at the Tevatron[[3] 
and the Next Linear Collider (NLC) ||, the top quark spin is decomposed in the Helicity 
basis. Recently, Mahlon and Parke || proposed the Generic spin basis and found that 
the "Off-diagonal" basis, a special case of the Generic spin basis, is a more optimized 
decomposition of the top quark spins for e + e~ colliders because the contribution from 
like-spin pairs of top quarks vanishes to leading order in perturbation theory. 

The 0(a s ) QCD corrections to the cross sections of e~e + — *> tt process were calculated 
in refs.||. In this letter, we first repeated the calculation of 0(a s ) QCD corrections to 
this process. We then present, in the Generic spin basis, a general formalism of one-loop 
radiative corrections to the spin correlations in the top quark pair production at the NLC, 
and calculate the 0(a s ) QCD corrections and give the numerical results. 

In the SM, we consider the process 

e"e + tt (1) 

at the NLC with -y/i = (0.4, 1) TeV. The tree level Ve~e + vertex can be written as 

^Vee = ^(e V L P- + e v R P + ) (2) 

where P± = (1 ± 7s)/2, and the SM values for these coupling factors are e\ R = — 1 for 
V = 7, ef = (2 sin 2 9w — 1)/ (2 sin 6w cos 6w) and e R = sin 9w / cos 9w for V = Z, and the 
9w is the Weinberg angle. 
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The general Vtt (V = 7, Z) coupling can be written as 

(3) 



where ti (t^) is the momentum of the outgoing top (anti-top) quark. 

In the ti center of mass frame (CMS), the scattering plane is defined to be the X-Z 
plane where the electron is moving along the +Z direction and 9 t is the scattering angle 
of top quark, and we also set 4>t = 0. The Born helicity amplitudes for the process ([!]) 
are obtained by summing the contributions from both the Z and 7: 

2 e 2 1 
M (h e -,h e+ ,ht, h) = -= M (h e -,h e+ ,h t , h) 1 + M (h e -,h e+ ,ht, hi) z R{s) . (4) 

V s 

where s = AE 2 is the total energy in CMS, E is the energy of the electron beam, and 
R{s) = s/(s - Mf ) Q. 

In the Generic spin basis, as illustrated by Fig.l in ref.0, the top quark (anti-top 
quark) spin states are defined in the top quark (anti-top quark) rest-frame, where one 
decomposes the top (anti-top) spin along the direction s t (sj), which makes an angle £ 
with the anti-top (top) momentum in the clockwise direction. Thus, the state tjtf (^j.) 
refers to a top with spin in the +s t (— St) direction in the top rest-frame, and an anti-top 
with spin +% (— % ) in the anti-top rest-frame. 

In the Generic spin basis, the amplitudes M (h e -, h e +, s t , Sf) for the process e~e + — > ti 
can be generally written as 

M ( — h tftf or t^j) = =f [m t A L sm9cos£ - (EA L cos9 + KB L ) sin£] , (5) 
M Q (- + t^i l or tfa) = -(EA L cos9 + KB L )cos^ -m t A L sm9sm^ 

T(EA L + KB L cos6), (6) 

M$(+ -t^ or t\ii) = =F [m t A R sm9cos£, - (EA R cos6 - KB R ) sin£] , (7) 
M (-\ — tft| or t^) = — (EA R cos9 — KB R ) cos£ — m t A R sm9 sin^ 

±(EA R -KB R cos6), (8) 



1 At NLC with y/s > 400GeV, the imaginary part of the Z propagator can be neglected safely. 
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where K = (E 2 — m 2 ) 1 ^ 2 . The amplitudes in Helicity basis can be obtained easily by 
setting cos£ = ±1 in Eqs.([5])-(||). The form factors Al,r and Bl R are defined as 

Al,r = ^(el R A y + el R A z R(s)), (9) 
Bl,r = ^(el R B, + el R B z R(s)). (10) 
At tree level, the form factors (Ay, By) (V = 7, Z) appeared in Eqs.(0),(H) an d fllCf) are 
A° = l B» = 0, A%= *: % r H Z , ^=_1__. (11) 



When we make a special choice for the angle £ in the Generic spin basis, 

m t A° r sin , , 

ta ^ = ^cosWj«r (12) 

the Generic spin basis turns into the so-called Off-diagonal basis || for e^e^ 
scattering 0, then we will have 

+ ort t) = (13) 

a cos 

d<7{) ' +t T t i0 r^ T ) = (-^_)|M (- + t T t i ort i t T )| 2 (14) 

+ -t^ort l t l ) « (15) 



dcos# v ! ! v 32tts 



(icos^ 

c/ctq / r Z \ f ^ ,1 \ r / ,7 7 1 1 2 



f/f , )S 0( + -*T*ior*i*T) = (^)l^o(+-^|Ort i t T )| i (16) 

where the subscript means the tree level physical quantities. 

At tree level, the form factors Ay and By (V = 7, Z) are very simple, and Cy = 
Dy = 0. When we consider the one-loop corrections, the form factors may become much 
more complicated. In general, including the one-loop corrections, the form factors can be 
written as 

Ay = A° v + 5Ay, B V = B V + 5By (17) 

where 5 Ay and 5By represent the one-loop corrections. 

2 As pointed in ref.0], there are two Off-diagonal basis for &fe\ and e^ej scattering respectively. But 
these two Off- diagonal bases are almost identical Q and we also only use the Off-diagonal basis for e.~fe^ 
defined by Eq.@ even when discussing the case of e. R e\ scattering. 
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At one-loop level, in the Generic spin basis, the general differential cross sections for 
e~e + scattering are 



(h e -,h e +,s h st) = tJ^— \M (h e -,h e +,s t ,St) + 5M(h e - , h e +, s t , s?)| 2 



dcos9 y c ! ' 32ns 
where Mq is the tree level amplitudes given in Eqs 
one-loop contributions, which are given by 



and the 5M represents the 



SM^+t^ort^) 



SM(+-t^ort l i l ) 



SM(+-t ] i l ort l i{ / 



=F (m t SA L - K 2 C L ) sin6»cos£ 

- (K5B L + E5A L cos 9) sin f ] - EKD L sin 0, (19) 

- (m t 5A L -K 2 C L ) sin 9 sin f - (K6B L + E5A L cos 0) cos £ 
=F(£<L4 L + A:55 L cos0) (20) 
=F (m t 8A R - K 2 C R ) sin cos £ 

-(£<L4 fl cos0 - #6-8^) sin£] - EKD R sin9 (21) 

- (m t <L4 B - K 2 C R ) sin 6 1 sin f - (-E^Ar cos 6 1 - ^55^) cos £ 

=F (K5B R cos - £5 (22) 



with 



5A L)R = 


Ae 2 E 


s 


8B LjR = 


Ae 2 E 


s 


Cl,r = 


4e 2 E 


s 




Ae 2 E 


D L ,R = 


s 



(el tR 5A^ + ef tR 5A z R(s 
(el R 5B, + el R 5B z R(s)) , 
(el, R C, + el R C z R(s)) , 
(elnD-y + el R D z R(s 



(23) 
(24) 
(25) 
(26) 

The above formulae are valid for the general cases corresponding to the vertex couplings 
as given in Eq.([3]). 

In what follows, we will calculate the 0(a s ) QCD corrections to the top quark pair 
production as an interesting example. The relevant Feynman diagrams include the QCD 
virtual corrections as well as the real bremsstrahlung graphs. 

We will use dimensional regularization to regulate all the ultraviolet divergences in the 
virtual one-loop corrections. To regulate the infrared divergences, we introduce a gluon 



mass parameter A in the gluon propagator, which is justified since the non-Abelian nature 
of QCD does not enter at the order in a s . We also adopt the on-shell mass renormalization 
scheme. 

For 0(a s ) QCD virtual corrections, the form factor Dy (V = 7, Z) is still zero, while 
the non-zero factors SAy, SBy and Cy are 

6A, = £-A°U(J3), 5A z = ^A z U(f3), (27) 

SB, = 0, 5B z = ^B%\U{(3)-2V{(3)i (28) 



C i = -^T C >^)i °z = -^ZZT^ Z V(P), (29) 



3irm t 7 " 37rmi 



with 



U{p) 



1 + /3 2 



+ 



3/?ln[X t 



ln[X t ] + 2 



M^H, (so) 
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where (3 = ^1 - 4m t 2 /s, ^(/5) = (1 - /5 2 ) ln[X t ]//5, X t = (1 - + /3), and Sp{z) = 

— JqIti (1 — t)/tdt is the Spence function. Since we have neglected the width of Z gauge 
boson, only the real parts of 8 Ay, 5 By and Cy will contribute to the corrections, and 
therefore we show only the real parts of them in Eqs.(|27|)-(|29|). In fact, the contributions 
from the imaginary parts of 5 Ay, 5 By and Cy are negligibly small even if we keep the 
imaginary parts of the Z boson propagator, by our numerical calculations. 

Although ultraviolet divergences have canceled in differential cross sections in Eq. fll8|) , 
the infrared- divergent part is still present. 

We regulate the infrared divergences associated with the soft real-gluon emission using 
the same gluon mass parameter A. Under the soft gluon approximation (SGA), the 
amplitude induced by the soft real gluon emissions can be written in a factorized form 
proportional to the tree level amplitude. The corresponding differential cross sections 
therefore are the form of 

d , aSOft Q (h e -,h e +, s t , St) = Vsoft J h e -,h e +,s t , %) (31) 
a cos 9 a cos 

The factor f] so ft can be written as 

Vsoft = Gir + Gs + Gf in (32) 



Here Gjr represents the infrared divergent part, G$ is the truncation part, and Gfi n is 
the finite part, and we have 



GlR 

Gg 

Gfi n 



3n s \ 



4a, , Auj, 



3tt 
4a, 



In 



HX t } + 



s3 

s — 2m\ 
s — 2m? 



2Sj9(l-X t ) + -ln 2 [X t ] 



(33) 
(34) 
(35) 



3n [P 1 sp 

where uo g is the fixed maximum gluon energy. 

All the infrared divergence will cancel after adding the real corrections in Eq. (|31~|) to 
the virtual corrections in Eq.fllgp, as they must. It is easy to see that the Off-diagonal 
basis defined at tree level is still valid after the inclusion of QCD virtual corrections and 
soft real gluon emission contributions. With the same value of £ as given in Eq. (|12D , the 
0(a s ) QCD corrected differential cross sections in the Off-diagonal basis are 



da 



dcos9 
da 

dcos6 



da 
dcos9 



da 
dcosO 
da 

dcos8 



- + trt 



da 



- + t ] t 



in, 



:+ - *t*t) 

;+ - hh) 



dcos9 

dao 
dcos9 

, KA L 
m t 

da 
dcos9 

, KA L 
m t 

da ^ 
dcos9 

dap { 
dcos9 



l L l) 







(36) 



0:.s 



3tt 



X{(3) + 2 



KV(/3) 



|M (- + t T t 1 



sin 9 sin £ — 2B L cos £ — 2B L cos 9) 



(37) 



3tt 



+ 2 



KV((3) 



|M (- + t i t Ty 



sin sin £ — 2fi^ cos £ + 2i?^ cos i 



:+-hf)li + ^- 



(38) 
(39) 



3tt 



X(/3) + 2 



KV{0) 



|M (+-t T t 4 



da 
dcos9 



(+ " ¥t) 



■( — sin^sin^ + 2_Br cos£ — 2_Br cosi 

c/cr 



(40) 



dcos9 

jKAr 

m t 



3tt 



X(/?) + 2 



W(/3) 



|M (+-t i f T 



sin sin £ + 25r cos £ + 2Br cos ( 



(41) 
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with 

"1 + /3 2 



+ 2 



(3 



H*t] + 2 



m 2 

In 



* + n Sp(1 _ Xt) _ i±^! hl|A ,| _ i+r lAXt] _ 4 , (42) 



/? L tJ p 

From above analytical expressions, we can see that the one-loop QCD corrections do not 
change the spin configuration of top quark pairs. For e^e\ scattering, the contribution 
from the up- up (UU) and down-down (DD) top quark pairs vanish at the 0(a s ) order in 
Off-diagonal spin basis. For the e^ej scattering, the contributions from the UU and DD 
pairs of top quarks are also very close to zero. The cross sections in Helicity basis can be 
obtained easily by setting cos£ = 1 in the above formulae. 

In the numerical calculation, we use the following parameters as standard input 0, 

a s (M z ) = 0.118, m z = 9l.l87GeV, sin 2 6 W = 0.2315, 

a = — , m t = 175GeV, T z = 2A91GeV. (43) 
128 

For the running of ct s , we use yfs as the renormalization scale. For the infra-red cutoff u g 
we take u g = (y/s — 2m i )/5 f| under the soft gluon approximation. 

In the Off- diagonal spin basis, only one spin component is appreciably non-zero: the 
up-down (tfti) component for ej j e + scattering and the down-up (t|tf) component for e^e + 
scattering, as illustrated in Fig.l. In the Off-diagonal spin basis, we have 

da _ , - \ (0.23 ~ 0.70) pb, at tree level, 

'e L e+ ^ t ti ) = J (44) 



dcosd ' [ (0.33 ~ 0.89) pb, with QCD correction 

and 

da _ , j (0.09 ~ 0.33) pb, at tree level, 

{e R e + -> t^ T ) = I (45) 



dcosd [ (0.16 ~ 0.43) pb, with QCD correction, 

for —1 < cos$ < 1 and s/s = AOOGeV. Above results show that, the QCD corrections 
make the differential cross sections of dominant spin components larger by ~ 30% com- 
pared to the tree level ones for both e~[e + and e] i e + modes. For -^i = 1000 GeV, however, 
the increase of the corresponding differential cross sections induced by the inclusion of 
QCD corrections is less than 3% as shown in Fig.2. 



Fig. 3 shows the fractions of the total cross sections for the dominant spin components, 
defined in both the Helicity basis and the optimal Off-diagonal basis for e^e + scattering, 
as a function of the total energy y/s. The fractions were defined as 

Run = T » ( 46 ) 

0~UD + 0~DU 

Rlr = — ■ (47) 

°~t L t L + cr tL i R + cr tR t L + o tR t R 

In the Off-diagonal spin basis, the dominant spin component tft± makes up the 99.8% 
(96.4%) of the total cross section for y/s = 400GeV (ITeV) after including the 0(a s ) 
QCD corrections. In the Helicity basis, however, the fraction for the dominant spin 
component tj/tn is only 53% at y/s = 400GeV, as shown by the lower two curves in Fig.4. 
Although this fraction for tiin will increase to ~ 83% at y/s = ITeV, it is still less than 
the corresponding ratio Rt r t l = 96.4% at y/s = ITeV. For the e~ R e + scattering we have 
similar results. 

When defined in the optimal Off- diagonal basis, the spins of the t and i pairs produced 
from polarized efe + and e^e + scattering are uniquely determined, and this specific fea- 
ture will not be changed after including the 0(a s ) QCD corrections under the soft-gluon 
approximation. Because the contribution to spin correlations in the top quark pair pro- 
duction from the hard-gluon emissions is numerically very small, we here present only the 
results under the soft-gluon approximation. 

For the sake of cross check, we also present the numerical results in the ordinary 
Helicity basis. Fig.4 shows the differential cross sections for the process efe + — > t^^i and 
e i? e + — > t L,Rt in the Helicity basis with y/s = AOOGeV, and the corresponding differential 
cross sections are clearly well consistent with those shown in Fig. 8 of ref.||. 

To summarize, we have calculated the 0(a s ) QCD corrections to the spin correlations 
in the top quark pair production at NLC in the SM. We start from the general forms of 
Vti vertex (V — 7, Z), derive out the amplitudes in the Generic spin basis, and give the 
general formalism of including one-loop corrections to the differential cross section in top 
quark pair production at the NLC, and finally calculate the 0(a s ) QCD corrections. We 
found that: 
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(a) In Off-diagonal basis, the 0(a s ) QCD corrections to e~[e + ( e^e + ) 
scattering process increase the differential cross sections of the dominant spin 
component t\t^ (t]t{) by ~ 30% and ~ (0.1%— 3%) depending on the scattering 
angle for ^/s = AOOGeV and ITeV, respectively. 

(b) The Off-diagonal basis remains the optimal spin basis even after the 
inclusion of 0(a s ) QCD corrections. At ^J~s = AOOGeV, the dominant spin 
components in both ej j e + and e^e + scattering make up more than 99% of the 
total cross section at both tree and one-loop level, but such fraction is only 
~ 53% in the Helicity basis. 

Note added. While preparing this manuscript the paper of J. Kodaira, T. Nasuno 
and S. Parke appeared || where the QCD corrections to spin correlations in top quark 
production at e + e~ colliders is also calculated. Their numerical results under the soft 
gluon approximation are in very good agreement with ours. 
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Figure Captions 

Fig.l: The differential cross sections in the Off- diagonal basis for the e~[ R e + — > tt pro- 
cesses UD (ffij, DU (tj,t T ) and UU + DD (i T i T + assuming y/s = AOOGeV 
and c<j 9 = (y/s — 2m t )/5 = lOGeV. The non-dominant (i-j-fj.) component for the 
e^e + (e^e + ) scattering is multiplied by a factor of 100. 

Fig.2: The same as Fig.l but for ^1 = ITeV and u g = (y/s - 2m t )/5 = ISOGeV. The 
non-dominant tjtj- + component for the e^e + scattering is multiplied by a factor 
of 10. 

Fig. 3: The y/s dependence of fractions of total cross sections for dominant t^t^ (t L t R ) 
component for ej j e + scattering. The two upper curves show the fractions in the 
Off-diagonal basis, while the two lower curves correspond to the fractions in the 
Helicity basis. 

Fig. 4: The differential cross sections in the Helicity basis for the e~[e + — > t LjR t and 
e~ R e + — > t LiR t processes at y/s = AOOGeV, and assuming u g = (y/s — 2m t )/5 = 
WGeV. 
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